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r <3-> ' ABSTRACT 



The broad absorption lines (BALs) of QSOs indicate abundances of heavy elements, 
relative to hydrogen, that are 1 to 2 orders of magnitude higher than the solar values. 
In at least one QSO, an especially large enhancement of phosphorus is observed. These 
abundances resemble those in Galactic novae, and this suggests that novae may produce 



the BAL gas. The needed rate of nova outbursts may come from single white dwarfs 



CN| ■ that accrete gas as they pass through a supermassive accretion disk around a central 



black hole. 
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BAL QSOs show broad absorption troughs attributed to rapidly outflowing gas, of unknown 
origin, located outside the continuum source and the broad emission-line region. BALs typically 
occur in La, C IV A1549, Si IV A1400, N V A1240, O VI A1034, and sometimes Mg II A2798 and Al 
III A1857 (see reviews by Weymann, Turnshek, and Christiansen 1985, "WTC"; and Turnshek 1988, 
1995). BALs often set in at the systemic velocity of the QSO, but in some objects "detached" BALs 
begin at a velocities up to ~ 10 4 km s _1 . The absorption often extends to velocities ~ 3 x 10 4 
km s , with varying degrees of structure and residual intensity through the line profile. The 
derived column densities for the absorbing ions (e.g., C +3 ) are ~ 10 16 cm~ 2 . BALs occur in about 
10 percent of radio quiet QSOs (Weymann et al. 1991) and rarely in radio loud QSOs . Either 
most radio quiet QSOs have BAL material covering about 10 percent of the sky as seen from the 
continuum source, or a subset of QSOs have BAL regions with a larger covering factor. The BAL 
region may have a disk-like geometry (Turnshek 1995; Goodrich and Miller 1995). The absorbing 
material must have a transverse dimension ;> 10 16 cm, sufficient to cover the continuum source. 
The N V BAL often strongly absorbs the La emission line, implying a radius ^ 10 18 cm (Turnshek 
et al. 1996). The absorbing clouds, which occupy only a tiny fraction of the volume, may be 
accelerated by radiation pressure involving the resonance lines or by ram pressure of a fast moving 
wind (WTC). Proposed sources for the BAL gas include winds from red giant stars (Scoville and 
Norman 1995) or from a supermassive accretion disk (Murray et al. 1995). 
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2. Chemical Abundances 

The BAL gas has remarkably high abundances of heavy elements, relative to hydrogen. The 
abundances are derived from column densities of ions (in turn derived from observed BAL optical 
depths) together with photoionization models of the ionization equilibrium. Junkkarinen et al. 
(1987) argued for a minimum (Si/H) > 10(Si/H)Q for their sample of BAL QSOs, and their discus- 
sion suggests a likely value (Si/H) ~ 25(Si/H) Q . Turnshek et al. (1987) reported a C/H ratio 10 to 
100 times the solar value in Q0932+501 (see also WTC); and for Q1413+113, they found that S/C 
and either P/C or Fe/C are at least 100 times solar. For Q0226-1024, Korista et al. (1992, 1995) 
found a heavy element abundance Z/Z rs 5 to 10, in the context of an assumed chemical evolution 
scenario. Analyzing the same observations of Q0226-1024, Turnshek et al. (1996) developed atwo 
component photoionization model that gave ratios ~ 4, 3, 9, and 120 times solar for N/O, N/Si, 
N/C, and N/H respectively (best fit model). Enhancements of other elements, including S and Ar, 
are observed in some cases (Turnshek 1995). Remarkably, P/C is ~ 65 times the solar value in PG 
0946+301 (Junkkarinen et al. 1995; Junkkarinen 1995). The heavy element enhancements can be 
reduced, but not eliminated, by using complicated shapes for the ionizing continuum (WTC). 

Novae show very large enhancements of C, N, O and sometimes Ne, Mg, Al, Si, S, Ar, Ca, and 
Fe (Andrea et al. 1994, and references therein). Novae are produced by thermonuclear explosions on 
the surface of white dwarf stars accreting hydrogen rich gas from a close companion star (see review 
by Starrfield 1989). The ejected gas is enriched with heavy elements mixed in from the white dwarf. 
Enhancements of Ne and heavier elements, observed in "neon novae", are attributed to events on 
O-Ne-Mg white dwarfs (Politano et al. 1995, and references therein). Hydrogen burning leads 
to especially large enhancements of nitrogen, and in the case of O-Ne-Mg white dwarfs, heavier 
odd numbered elements. Fig. 1 illustrates the abundances in a set of novae, determined in a 
uniform manner by Andrea et al. (1994). Abundances of C, O, and Si are typically enhanced by 
one to two orders-of-magnitude, compared with solar values; and N/H is enhanced by two to three 
orders-of-magnitude. 

I have estimated abundances for several BAL QSOs by assuming that their elemental abun- 
dance ratios scale from those given by Turnshek et al. (1996) for Q0226-1024 in proportion to the 
corresponding ionic column densities. (This assumes that the ionization corrections are similar 
in all these objects. Column density ratios for ions of a single element are generally unavailable 
to constrain models for individual objects.) Column densities were taken from Junkkarinen et al. 
(1987) for 5 BAL QSOs and from Junkkarinen (1995) and co-workers for Q0946+301. Figure 1 
shows the resulting abundance ratios along with a value (C/H) = 30(C/H) Q for Q0932+501 (Turn- 
shek et al. 1987). The abundance distributions for C, N, O, and Si for novae and BAL QSOs are 
similar. Politano et al. (1995) predict P/C ratios of ~ 50 and 300 times solar for nova explosions 
on O-Ne-Mg white dwarfs with masses of 1.25 and 1.35 M , respectively, bracketing the value for 
PG0946+301. These results suggest that the BAL gas may come from novae in a massive star 
cluster in the active galactic nucleus. 
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High Si/C and P/C in BAL QSOs suggests that neon novae are typically involved. If O-Ne-Mg 
white dwarfs come only from progenitors of mass 8 — 10 M (Nomoto 1984), they should be only 
a few percent of Galactic white dwarfs. Truran and Livio (1986) attempted to explain the high 
observed incidence of neon novae in the Galaxy in terms of a small accreted mass per outburst. 
However, observed masses of neon nova shells are ~ 10~ 4 M (e.g., Shore et al. 1993), so another 
explanation may be needed. If O-Ne-Mg white dwarfs come from progenitors down to ~ 5 M Q 
(Chiosi et al. 1989), these would be a large fraction of all white dwarfs for a cluster age ~ 10 8 yr, 
the likely duration of a QSO episode (Norman and Scoville 1988). Alternatively, perhaps nuclear 
burning to high atomic numbers, which Politano et al. (1995) find only for the most massive white 
dwarfs, actually occurs for lower masses as well. Finally, if the BAL geometry is such that debris of 
different novae are comingled, as might occur in the disk geometry of Section 4, neon novae could 
contribute heavier elements to the mix. 

The heavy element abundances in the broad emission-line gas of QSOs may be up to an order 
of magnitude higher than solar, relative to hydrogen (Hamann and Ferland 1993), but apparently 
are not so high as in the BAL gas. The emission-line abundances have been attributed to rapid 
chemical evolution involving the usual stellar sources of heavy elements (Hamann and Ferland 
1993), and this would not explain high P/C. This suggests that the two regions likely have different 
origins. 



3. Ordinary Novae 

Novae are minor sources of interstellar gas in the Galaxy. Can they nevertheless be the dom- 
inant source of BAL gas? We first consider ordinary novae in a massive nuclear star cluster and 
then the possibility of accretion onto single white dwarfs passing through a supermassive accretion 
disk. 

The mass of a nova shell typically is M e j « 10~ 4 3 M (Warner 1989). As the shell expands and 
accelerates radially outward, it presumably fragments into a collection of clouds or filaments. We 
assume that the lateral expansion continues at the original ejection velocity, v nov « 2000 km s _1 , 
appropriate for neon novae (e.g., Ferland, Lambert, and Woodman 1977; Shore et al. 1993; Gehrz 
et al. 1985). The column density of carbon atoms through the debris is Nq ~ Mq / (irr^ h m c ) , where 
r s h is the shell radius, Mq is the mass of carbon in the nova shell, and mp is the mass of a carbon 
atom. Therefore, we may write 

r sh « (10 17 - 3 cm)My%N-%\ (1) 

where M c ,-5 = M c /W~ 5 M and 7V"c,i6 = N c /W 16 cm . For a carbon abundance C/H = 
10 1 - 5 (C/H) Q , we expect M c ~ 10~ 55 M . Assuming that roughly half of the carbon is C +3 , we 
find Nc,i6 ~ 1 f° r r sh ~ 10 170 cm. This radius is large enough to cover the continuum source and, 
within the uncertainties, the La emitting region. As the debris expand, the continuum radiation 
pressure presumably ablates off material in the form of small clouds, which move radially outward 
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as they are accelerated to the observed outflow velocities (cf. WTC). For AR/r s h ~ w/v nov , Nc,w 
drops to unity for distances A.R w (w/v nov )r s h ~ 10 18 cm from an initial location at radius Ri n u, 
where 10, 000 km s 1 is the outflow velocity. 

The expected covering factor for nova shells depends on the nova rate in the star cluster in 
the galactic nucleus. The covering factor for one shell is £l s h/4n ~ 7rr 2 ft /(47ri? 2 ) ~ 0.25(v nov /w) 2 , 
where the second equality assumes Ri n u < AR. The number of shells in play at a given moment 
is N s h Ri NR/w, where N is the rate of nova explosions and R/w ~ r s h/v nov is the crossing time. 
We evaluate Q/4tt for the value of r s h that gives Nc,ie ~ 1. This leads to 

fl/47r « 0.25Nr sh v nov w- 2 « lO " 15 ^^/™)^ 1 Mc%Nc$, (2) 

where iV is the nova rate per year and = w/(10 9 cm s _1 ). If Rinit Ai£, this result should 
be adjusted accordingly. Note that Q/4tt varies as (N 'M e j)M~^ 2 , so that events involving a small 
mass, such as novae, are highly effective for a given average mass- loss rate, NM e j. The nova rate 
in the Galaxy (mass ~ 10 11 M Q ) is ~ 40 yr~ 4 (Warner 1989, Ciardullo et al. 1990); and therefore 
we take N a 10~ M c / 8 , where M c / ]8 is the mass of the nuclear star cluster in 10 8 M . Then we 
find 

0/47T ~ lQ~ 1A M 8 (v nov /w)wg 1 » 10" L8 M 8 . (3) 

A ten percent covering factor requires Ms ~ 10. Analogous estimates of the covering factor for 
BAL absorption due to planetary nebulae, supernova remnants, and "stellar contrails" (Scoville and 
Norman 1995) indicate that novae are competitive or dominant by up to an order of magnitude. 
However, none of these sources match the BAL abundances as naturally as do novae. 

The mass loss rate in carbon alone can be estimated as 

M c ~ N c m c ^R 2 (R/wy 1 (n/A7:) w (10" 4 - 4 M Q yr^Nc^Risw^Q/ ^). (4) 

For iVc,i6 ~ #18 ~ ^9 ~ 1 and 17/4vr 0.1, this gives M c ~ 10~ 5 - 4 M yr^ 1 . For C/H 
30{C / H)q, the total mass loss rate then is ~ 10~ 4 ' 2 M Q yr _1 , a rather modest value. 

Norman and Scoville (1988) discussed a coeval star cluster of mass ~ 10 9 Mq whose evolu- 
tionary debris fuel the central black hole. The age of the nuclear star cluster is at most 10 9 ' 3 yr 
at redshift z = 2 for q Q = 1/2 and H~ l = 15 Gyr, and it could be as young as the estimated 
QSO lifetime of ~ 10 s yr. The "evolutionary flux" of stars leaving the main sequence drops an 
order-of-magnitude as a coeval cluster evolves from age 10 8 to age 10 9 ' 5 years and somewhat fur- 
ther by age 10 10 years (Norman and Scoville 1988). The effect of this on the nova rate is unclear. 
The required star cluster, confined to a radius not much larger than 10 18 cm, would have a stellar 
velocity dispersion of several thousand km s~ 4 . The stellar collision time would be shorter than 
the Hubble time and possibly even the estimated ~ 10 8 yr lifetime of a QSO episode. Observations 
of nearby galactic nuclei offer little support for such massive nuclear star clusters (e.g. Lauer et 
al. 1992, and references therein). Thus, we are motivated to consider ways to enhance the nova 
rate in QSOs. 
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4. Single White Dwarfs 

An intriguing possibility is suggested by the work of Artymowicz, Lin, and Wampler (1993, 
ALW) . They consider stars on orbits passing through an accretion disk around a central black hole 
of mass M/j, accreting disk gas during each passage. White dwarfs orbiting through the disk will 
also accrete disk material, and this raises raises the possibility of nova explosions on single white 
dwarfs that have accreted the requisite amount of hydrogen rich gas (cf. Truran et al. 1977). The 
circular velocity in units 10 8 km s _1 is v c ^g ~ 10 ' 1 M^gR^^ 2 , where Mtot = + M c \. The 
average stellar velocity, u*, will be roughly v c . Let v re i be the velocity of a star relative to the 
orbiting material in the disk. If the cluster is corotating with the disk, most stars will have v re \ 
substantially less than v*. For parameters of interest, v re i is less than the escape velocity from the 
surface of a white dwarf and large compared with the sound speed in the disk. The accretion rate 
onto the white dwarf while in the disk is 

m pa 2.57rG' 2 m 2 /w~J (5) 

(Bondi and Hoyle 1944), where m ~ 1 Mq is the mass of the white dwarf and p is the gas density 
in the disk. Following ALW, we assume that the disk thickness H is such that the Toomre (1964) 
stability parameter, Q, is near unity. Consequently, itY,R 2 /Mh ~ H/R pa 1CT 2 , where £ = 2Hp 
is the disk surface mass density and a value H/R ~ 10 -2 corresponds to the disk's expected 
vertical equilibrium. Then £ pa (10 2 - 8 g cm^M^isRrg 2 , P ~ (10~ 13 ' 5 g cm _3 )M/ li8 i?j~ 8 3 , and m pa 
10~ 10A Mh t 8Ri8 v rei 8' w h ere m = rh/1 M yr _1 . If the star's vertical velocity v z is ~ v re i/\/3, 
then the duration of passage through the disk is At pa (10 7 yx)Ri8V~^ 8 , and the mass accreted 
is Am pa (10~ 9,4 M^Mh^Ri^v^ g . If this much mass is accreted twice per orbital period, P m 
iTxRjv^ = (10 3 - 3 yr)Ri8V~l, then the average accretion rate is (m ) m (10^ 12A )Mh t 8Ri8 v *,8V~^i 8 - 
A typical white dwarf will undergo an explosion at intervals t nov pa (10 8 - 1 Y^)M^R\ % v~^v^. el g . 

Suppose there are N\vd ~ 10 7 M c / i s white dwarfs in the nuclear star cluster (Allen 1973). (For 
a Salpeter [1955] initial mass function ranging from 0.2 to 20 Mq and a main sequence turnoff of 
1.4 Mq, corresponding to an age of 10 9 ' 3 yr, one has Nwd pa 10 6 ' 8 M c / §.) Let the stars have a typical 
velocity a < in the frame corotating with the disk. Then the nova rate for most white dwarfs, 
with v re i ~ cr, will be N Q Pa lO-^M^sMd^R^v^s 4 - With M h ~ M d and v 2 pa GM tot /R, this 
gives N pa 10~ 1 - 9 M t 1 o { 2 g i? ] ~ 8 3 ' /2 (^/(T) 4 . For i?i 8 pa M 8 pa 1 and cr/w* pa 0.5, we then have 

Q/iir pa 10-°- 8 (v nov /w)w 9 1 pa 10" 1 ' 5 . (6) 

However, accretion is most rapid onto stars with small v re i, and these can dominate the nova rate 
in spite of their relatively small numbers, N(v re i)/N pa (v re i/a) 3 . If relative velocities down to 
some minimum, v m i n , contribute to the nova rate, then a simple integration over v re i shows that 
the preceeding expressions for iV and 0/47T should be increased by a factor ~ (cr/v m i n ). If violent 
nova explosions occur for m < 10 -7 M Q yr -1 (van den Heuvel et al. 1992), then we may take 
Vmin ~ 10 6 ' 6 cm s _1 pa lO^cr. This increases N Q and Q/4tt by an order of magnitude, so that 
S7/47T pa 10 -0 ' 5 . Within the uncertainties, this is consistent with the observed 10 percent incidence 
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of BALs. For the assumed parameters, novae involving single white dwarfs outnumber ordinary 
novae by a factor ~ 40. 

Indications of a flattened BAL region are consistent with novae from white dwarfs with small 
v re i/v*, which would orbit close to the disk (z/R ~ v re i/v c ). The actual nova event would occur at 
a random place in the white dwarf's orbit. The association of BALs with radio quiet QSOs might 
involve the absence of a suitable accretion disk in radio loud objects, perhaps because the latter 
typically occur in elliptical galaxies. 

5. Discussion 

In summary, the high chemical abundances derived for the BAL region of QSOs resemble the 
abundances of Galactic novae. The mass of heavy elements in a nova shell is consistent with con- 
straints on the dimensions and location of the BAL absorbing material. A very massive nuclear star 
cluster must be postulated in order for ordinary nova outbursts to explain the observed incidence of 
BALs in QSO spectra. Accretion of gas by white dwarfs passing through a supermassive accretion 
disk can provide sufficient nova outbursts with a more modest star cluster. 

Observational consequences include a small BAL region that could show changes over a few 
years. Different velocity components in the BALs of a given object could result from different novae 
and show different abundances. Odd numbered elements such as Al and P should be especially 
enhanced in abundance. 

This work was inspired in part by the talks by P. Artymowicz and V. Junkkarinen (who 
mentioned high P in nova models) at the June, 1995, meeting of the American Astronomical 
Society. I am grateful to F. Bash, R. Ciardullo, D. Lin, K. Korista, D. Lambert, S. Starrfield, D. 
Turnshek, and J. C. Wheeler for valuable discussions, and to V. J. and D. T. for communicating 
results in advance of publication. This work was supported in part by STScI grant GO-5434.03-93A. 
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Fig. 1. — Histogram of C, N, O, and Si abundances relative to H for novae and BAL QSOs, in bins 
of 0.5 dex, normalized to solar ratios. In each panel, novae are shown above the horizontal axis and 
BAL QSOs below it. Each small tick on the vertical axis represents one object. Results for novae 
are from Andrea et al. (1994), and those for BAL QSOs are from several sources described in the 
text. 



